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ABSTRACT 

There has been mounting observational evidence in favour of Luminous Blue Variables (LBVs) being the direct progenitors of su- 
pernovae. Here we present possibly the most convincing evidence yet for such progenitors. We find multiple absorption component 
P-Cygni profiles of hydrogen and helium in the spectrum of SN 2005gj, which we interpret as being an imprint of the progenitors 
mass-loss history. Such profiles have previously only been detected in Luminous Blue Variables. This striking resemblance of the 
profiles, along with wind velocities and periods consistent with LBVs leads us to connect SN 2005gj to an LBV progenitor. 
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In the canonical pic ture of the evolut ion of very massive (M > 
40 M Q ) stars (e.g. lLanger et al.lll994l) . the role of mass loss is 
widely accepted as being the key process which drives O and 
B-type stars into the luminous blue variable (LBV) phase, be- 
fore entering the He-burning Wolf-R ayet (WR) stage, at the 
end of which, lasting a t least 10 5 yrs dMevnet & Maedet1 l2003; 
Eldridge & Vink 2006), the object is expected to explode as a 
core-collapse supernova (CcSN). 

Of late, there has been a growing body of observational ev- 
idence that points to LBVs as the direct p rogenitors of a subset 
of CcSNe. Recently iKotak & Vink) (120061) suggested that LBVs 
themselves might explode; they proposed that the quasi-periodic 
modulations observed in the radio light curves of some SNe were 
a manifestation of variable mass loss during the S Doradus phase 
of LBVs, during which the star exhibits brightness variations 
of ~l- 2 magnitudes. This variability i s commonplace amongst 
LBVs (Hum phreys & Davidson! fl 994) and is distinct from the 
giant outbursts such as those of P Cygni or 77 Carina. 

Interestingly, the progenitor of SN 2006jc was observed 
to have a gia nt eru ption just two years prior to explosio n 
dNakano et id] 120061: Foley et al.l 120071: iPastorello et aD l2Q07h . 
Fole y et alj d2007l) argue that the progenitor star was of an early 
type (WNE) Wolf-Rayet star. Given that WR stars have never 
previously been obse rved to undergo an LBV-like eruption, 
Pastorello et al.l (120071) suggest a massive binary system consist- 
ing of an LBV and a WR star as an alternative progenitor system 
for SN 2006jc. In this scenario, the pre-SN explosion coincident 
with SN 2006jc is attributed to the LBV, while the SN resulted 
from the WR companion. A simpler explanation would be to ac- 
cept that a single, presumably massive star exploded during or at 
the very end of the LBV phase. 
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There are other indications that LBVs may explode. 
iGal -Yam et al. I (120071) detected a very luminous source in pre- 
explosion images of SN 2005gl, but it remains to be confirmed 
whether or not this is a single object. Other interesting hints 
come from the similarities in LBV ne bula morphol ogies and 
the circumstellar medium of SN 1987A dSmith Il2007ab . Finally, 
the luminous type Iln SN 2006gy, ma y also have undergo ne an 
77 Car-type eruption prior to explosion ( Smith et a l. 2007b). 

A potential problem with a scenario in which LBVs ex- 
plode is that current stellar evolution models do not predict core- 
collapse during or soon after the LBV phase; indeed, in most 
evolutionary models, the star has not even reached its core He- 
burning phase. In other words, in current mod els, the core is 
not ev olved enough for core-collapse to occur. IPastorello et al.l 
(|2007|) point out that known LBVs that have undergone outbursts 
still retain some hydrogen and helium while SN 2006jc showed 
no spectral features due to hydrogen before ~65 d. However, it is 
not inconceivable that the pre-explosion outburst removed most 
of the outer layers, and/or that residual amounts of H or He may 
not necessarily give rise to strong lines. 

Here, we present the intriguing case of SN 2005gj for which 
we show spec t rosco pic evidence to support the suggestion of 
IKotak & Vinkl (120061) that the variable mass-loss of S Dor vari- 
ables betrays them as the endpoints of very massive stars. 



2. Data acquisition and reduction 

Optical spectra of SN 2005 gj were obtained on the 17 Dec. 2005 
and 1 Oct. 2006, with the Ultra- Violet Echelle Spectrograph 
(UVES) on the Very Large Telescope. These epochs correspond 
to 86 and 374 days respectively, as suming an explosion date of 
22 Sep. 2005 dAldering et alj|2006l) . Dichroic #2 was used to 
split the light into the blue and red arms of the spectrograph. 
The standard setup with central wavelengths at 4370 and 8600 
Awas employed, thus providing coverage in the blue of 3760- 
4985 A, and in the red arms of 6705-8520 A and 8660-10255 
A. For both epochs the target was observed for 2300 s with a slit 
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Fig. 1. Left: Optical spectra of SN 2005gj at 86 and 374 days after explosion, after correction for the redshift of the host galaxy. 
Right: Multiple component Gaussian profile fits to the Ha profile for both epochs; day 86 in top panel and day 374 in bottom panel 
(Table[TJ. Units of the axes are the same as the left panel. 



width of 1'.'6 and a seeing of ~ 1'.'05 providing resolutions of ~6 
and 4.5 km s in the blue and red arms, respectively. 

The data were reduced with the UVES pipeline implemented 
in the ESO-MIDAS software. Wavelength calibration was car- 
ried out using Th-Ar arcs taken on the same night as the 
SN 2005gj exposures, while flux calibration and telluric subtrac- 
tion of the data was done with respect to the standard LTT 1020. 
The flux calibrated spectra extracted from each chip and cor- 
rected to the rest frame are presented in Fig. [T| A redshift of 
0.0616 + 0.0002, as found bx Rldering et al.1 d2006l) . was adopted 
as it was consistent with our data. 

3. Spectral Analysis 

The optical spectra of SN 2005gj are dominated by Ha in emis- 
sion. The profile consists of a strong narrow emission compo- 
nent superimposed on broad emission wings (Fig. Q]). Such Ho- 
profiles are the hallmarks of SN Iln spectra ('n' standing for 
'narrow'), believed to arise from massive progenitors that have 
under gone substantial mass-loss prior to explosion (Schlegel 
1990). The multi-component profiles are characterised by broad 
underlying emission due to SN ejecta, an intermediate emission 
resulting from shocked material at the interaction front between 
the ejecta and the CSM, while the narrow component presum- 
ably arises from unshocked CSM photoionised by the SN. 

A striking feature of the P Cygni-like component of the Ha 
profile at 86 d is the presence of two clear absorption troughs 
(Fig. []]). This has previously not been observed in SN spectra. 
That this feature is associated with the SN is c lear, given the 
evolu tion of the profile between the two epochs. Alderin g et alj 
(2006) noted that the absorption component in the narrow P- 
Cygni profile was asymmetric and extended far out towards the 
blue edge. Due to the higher re solution of our spectra (c.f. resolu- 
tions of a few hundred km s _1 Aldering et al. 2006; Prieto et ail 
I20071) . we were able to resolve this narrow component. Careful 
inspection of the data revealed clear double absorption troughs 
in our day 86 spectrum associated with Hy, He 1 7065 A (Fig.|2]i, 
and also H6, and He i 6078 A albeit with poorer signal-to-noise 
ratios. By ~400d, the Ha profile had changed significantly: the 



Table 1. Parameters of the gaussian fits to the narrow compo- 
nents of H %a & He i 7065 A lines of SN 2005gj at 86 and 374 
days after explosion. 'Narrow E' - the narrow emission compo- 
nent, 'NAC (a)' and 'NAC (b)' - the blue and red absorption 
components respectively. Each profile has been fit by multiple 
gaussian profiles, representing the broad, intermediate and nar- 
row components of the profile. 



Fit properties 




Narrow E 


NAC (a) 


NAC (b) 


Day 86: H„ 


A (A) 




6564.29+0.02 


6558.35+0.06 


6561.25+0.06 


FWHM (km s~ 


') 


128 ± 2 


70+6 


133+21 


L (10 40 ergs s" 


) 


1.61±0.02 


0.14+0.01 


0.20 ±0.02 


Vedge (kms~') 






295 ± 57 


120 ± 75 


Day 86: H y 


A (A) 




4341.39+0.02 


4338.01+0.10 


4340.04±0.10 


FWHM (km s~ 


') 


97 ±3 


68 + 5 


115 ±4 


L(10 40 ergs s~ 


) 


0.15±0.04 


0.04+0.01 


0.12±0.05 


Vedge (km s _I ) 






286 ± 57 


115 ±64 


Day 86: Hei7065 


A (A) 




7066.33+0.03 


7062.61+0.09 


7063.97±0.10 


FWHM (km s~ 


') 


64 + 3 


22 ±9 


62±15 


L (10 40 ergs s~ 


) 


0.10+0.02 


0.006+0.001 


0.022±0.004 


Vedge (kms -1 ) 






160+ 59 


95 ±57 


Day 374: H tt 


-1(A) 




6562.60+0.04 


6559.70 ±0.12 




FWHM (km s~ 




88+ 10 


110 ± 23 




L(10 40 ergs s" 


) 


0.18+ 0.09 


0.028+ 0.005 




Vedge (km S~') 






190 ± 64 





intensity dropped by 17xl0~ 16 ergs -1 crrT 2 s , a single ab- 
sorption trough remained, and a broad blue shoulder had de- 
veloped, probably corresponding to the emergence of [O i]. To 
obtain quantitative information from these profiles we have de- 
composed them with multiple Gaussian profiles for the individ- 
ual emission and absorption peaks. The fits are shown in Figs.Q] 
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Fig. 2. Multiple Gaussian profile fits to the 86 d spectrum for Hy 
(left), and Hei 7065 A (right). The flux on the y-axes are given in 
units of 10~ 16 ergs -1 cm -2 s _1 . 
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& |2] and their physical parameters are summarised in Table Q] 
The Ha emission profiles at both epochs were fit by 3 Gaussians 
representing the aforementioned broad, intermediate and nar- 
row emission components, however this intermediate component 
was not detected in the Hy & He i 7065 A profiles. In addition 
to the emission peaks, the Ha profile from day 86 was fit with 
double absorption components, whilst only a single absorption 
trough was required for the day 374 spectra. 

Due to our incomplete wavelength coverage, particularly to 
the blue of Ha, the velocities determined from the Gaussian 
fits to the broad emission peak are subject to large uncertain- 
ties. From the B aimer lines we determine an ejecta velocity of 
30000 + 5 000 km s" 1 at day 86, which is consistent with that 
derived bv lAldering et alj from their day 1 1 spectra. At day 374 
the velocity decreased to 15500 + 5000 km s~'. The intermedi- 
ate Ha emission component was fit by a Gaussian with FWHM 
2560 + 30 km s -1 on day 86 and 1680 + 30 km s -1 on day 374. 
Whilst this measurement is a good indication of the shock veloc- 
ity for the symmetric profile presented at day 86, this is not the 
case for the asymmetric profile on day 374. In this later spectra 
the narrow Ha feature is redshifted with-respect-to the interme- 
diate component, this has been seen in a number of Type Iln 
spectra such as SN 1997ab, 1997eg & the hybrid Type Ia/IIn 
SN 2002ic dKotaketal .1120041) . This is generally attributed to 
self-absorption in the intermediate component which preferen- 
tially attenuates the red wing of the profile (Salamanc a et alj 
1998). A better measurement of the shock velocity is therefore 
the velocity of the blue edge of the profile at zero intensity, 
which is 2850 + 200 km s~' . The luminosities inferred from the 
na rrow line intens ities at day 86 agree well with those derived 
by lAldering et all with a significant decline in the luminosi- 
ties derived from the narrow component from 1.6xl0 40 ergs -1 
to 1.8xl0 39 ergs -1 by day 374. 

Other, less dramatic features in the spectra comprise a broad 
feature at ~ 8500 A, probably a blend of the Ca n infrared triplet 
with 1 8446. There is a weak, broad feature centered on 7300A 
which, depending on the SN type, can be identified as either [Ca 
n] or [Fe n], with the latter being typically present in the late- 
time spectra of SNe la. The [Ca n] lines however, are typically 
observed in the spectra of Type Iln and Ib/c SNe at >100d. 

4. SN 2005gj: Properties of the unshocked CSM 

These narrow P-Cygni features which are commonly found 
amongst Type Iln, are thought to represent the outflow of un- 
shocked CSM surrounding the SN event. As such, it provides 
insight into the environment prior to the SN explosion and on 
the wind properties of the SN progenitor itself. Stellar wind ter- 
minal velocities, Vco, are usually measured from the blue edge of 



Fig. 3. Comparison of the Ha profile of SN 2005gj in velocity 
space with that of the lu minous blue variables, A G Car & HD 
160529. Data taken from lStahl et all (1200 lll2003l) . 



strong resonance lines (v et jg e ). Unfortunately, there are no such 
resonance lines in the optical, and so we will use the Balmer and 
He i lines to estimate this velocity from their absorption com- 
ponents. As these lines are not sa turated, v e dgp is li kely to un- 
derestimate the terminal velocity dPrinia et alj [l990). From the 
Ha narrow absorption feature, we find a v e d g e of 295 km s -1 at 
86d with an additional component at 120 km s -1 whereas the 
component at 374d has an intermediate velocity of 190 km s . 
The Hy feature at day 86 gives velocities consistent with those 
found from Ha, however the bluer component of the He i profile 
appears at a lower velocity with a maximum v e d g e of 160 km s -1 . 

We can attempt to derive a n estimate of the mass-l oss rate 
from the Ha feature. Following ISalamanca et al.l d!998l) . L^ 1 = 
6h„MVj/4v w where, en a is an efficiency factor that peaks at ~0.1. 
Ljj' is the luminosity of the intermediate component, v w is the 
wind velocity and v s is the shocked CSM velocity. At day 86, 
this implies M = 6.4 and 2.6 xlO -2 M yr" 1 from the high (295 
km s -1 ) and low (120 km s -1 ) velocity components. At day 374, 
the derived M is 1.7x1 -2 M yr -1 using v w = 190 km s _1 and 
v s = 2850 km s , one should note that the shock velocity is a 
little less certain in this case (see Sect. [3]). These are very large 
mass-loss rates but are typical of those derived from the Ha lu- 
minosity in supernovae which show strong int eraction between 
the ejecta and CSM (e. g. SNe 199 7ab, 2002ic ISalamanca et all 
1998; Kot ak et al.ll2004l) . However. iKotak et al J found that mass- 
loss rates derived from dust measurements in the infrared are a 
magnitude lower. A similar observation can be made by com- 
pa ring the ma ss-loss rates of SN 2005gj derived here and that 
of iPrieto et alj from x-ray luminosities. Several caveats must be 
borne in mind when estimating mass loss rates via Ha profiles. 
Two of the most serious shortfalls involve the unknown ioni- 
sation conditions and uncertainties in e^ a . Realistic constraints 
can only be derived from detailed modelling which is beyond 
the scope of this work. 

5. Discussion 

SN 2005gj has been compare d to the prototype fo r the hybrid 
Type Ia/IIn objects, SN 2002ic ( lAldering et al.l2006l) . SN 2002ic 
was initially classified as a Type la, due to the presence of S 
n lines and a blue shifted absorption feature of Si n 6335 A, 
but was quickly reclassified as a hybrid Type Ia/IIn due to the 
presence of multi-component Balmer lines signifying interaction 
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with the surrounding circumstellar medium ( lHamuv et al.l2 003). 
This was an exciting discovery as the presence of circumstellar 
hydrogen in the spectra of a type la supernovae would shed light 
on the nature o f Type la progenitors. One scenario put forward 
bv lHamuv et al.l is that the dense CSM could represent mass-loss 
from an asymptotic-giant-branch (AGB) star. Unfortunately, the 
true classification and int erpretation of SN 2002ic remains un- 
der debate. In particular, iBenetti et al.l d2006l) argue that it can 
equally well be compared to a core-collapse SN, such as the 
type Ic SN 2004aw. This comparison is most convincing at late 
epochs, in fact the resembla nce of SN 2002ic to type Iln SNe 
had been noted earlier (e.g. Hamuy et al. 2 0031) . 

In spite of the Ha feature indicating CSM interaction, the 
early time optic al spectra of SN 2 02ic were undeniably Ia- 
like (see Fig. 3 lHamuv et al.l l2003l). lAldering et ail (120061) & 
iPrieto et al.1 (120071) argue that SN 2005 gj could be compared to 
the overluminous Type la, SN 1991T 'diluted' by the presence 
of significant CSM. However, the typical la spectral features due 
to S n and Si n which formed the crux of the la-evidence for 
SN 2002ic, are barely discernible in the spectra of SN 2005 gj. 
They attribute this to stronger interaction of the SN 2005gj with 
the CSM than for SN 2002ic, arguing that this claim is supported 
by its broader light curve and higher Ha luminosities. 

On the basis of the data presented here, we cannot verify 
or exclude the type la identification of SN 2005 gj. However, 
we present an alternative scenario that is consistent with the 
observations. We focus our attention on the double absorption 
troughs seen in the day 86 spectrum. Only in one other SN spec- 
tra, SN 1998S, hav e double absorption troughs been dectected 
dBowen et al.l 12000; Fas sia et aT]|200l|) . However the shape of 
the lines differ significantly. In SN 1998S a narrow P-Cygni pro- 
file of ~50 km s is superimpos ed on a broad, shallow absorp- 
tion profile dBowen et al.l l2000l see Fig. 8). Additionally, this 
absorption decreases in relative intensity from lower to higher 
Balmer lines, suggestive of an optical d epth dependence tha t is 
inconsistent with the SN 2005 gj spectra. IChugai et all d2002l ) in- 
voke several mechanisms that could accelerate the circumstellar 
gas around SN 1998S. These models have been tailored to match 
the early-time spectra of SN 1998S and it is difficult to see how 
these might be tailored to the significantly later epochs of our 
data. This type of test may provide some insight but is beyond 
the scope of this letter. 

Here we would like to stress the likeness of the profiles in 
SN 2005gj with another stellar object: Luminous Blue Variables. 
Multiple-absorption component Ha profiles have been detected 
in many LB Vs, viz AG Carinae, R127, R81, H P 1 60529 R66 
& P-Cvgni dWolf et alJfT98Tt IStahl et all fl983llaL l200ll 12001 
iLeithereret al.lll994l) . Fig. 3 shows several examples of the Ha 
profiles of LBVs gleaned from the literature. In LBVs these 
multiple-absorption compo nent profiles sign i fy the ir history of 
episodic mass-loss events. IVink & de Koterl d2002l) interpreted 
the variable mass-loss and wind velocity behaviour of LBVs as 
being due to a change in the ionization of the dominant wind 
driving ion Fe (wind bistability). This mechanism was also in- 
voked to explain the variations in the radio lightcurves of tran- 
sitional SNe (Kotak & Vink 2006). However, regardless of the 
mechanism the resemblance of the double-troughed LBV Ha 
profiles to that of SN 2005gj remains. 

An inter e sting hydrodynamical simulation by 
Ivan Marie et al.l d2007l) shows the evolution of CSM around a 
6OM star as it evolves from main-sequence to core-collapse via 
an LBV and WR phase. This simulation shows the formation 
of multiple absorption components from CSM shells formed 
by a fast WR wind running into an LBV shell. However these 



multiple shells occur in the very early stage of the WR-lifetime 
and hence if such a scenario was implied for SN 2005 gj, its 
WR-phase would have been short-lived. 

While the morphological resemblance of the line profiles 
observed in SN 2005gj and in LBVs is impressive enough, 
an even more compelling similarity is that of LBV wind ve- 
locities with those derived from the absorption components of 
SN 2005gj. LBV velocities lie in the -50-300 km strange, in 
strikingly good agreement with those derived from SN 2005gj 
(Table 1), which are larger than those observed in AGB winds 
(~10kms _1 ) and red-supergiant winds and much lower than 
typical WR-star winds. Typical timescales of S Doradus mass- 
loss events in LBVs are on the scale of years up to ~100 
yrs dvan Gendered l200lh . From a simple calculation we can 
determine the timescales of the profile variations detected in 
SN 2005gj. Using an intermediate ejecta velocity of 22500 
km s , a time variation of 288 days and an edge velocity of 300 
km s 1 we get P w =R e jecta/v e dge ~ 60 yrs. Thus the timescales 
associated with mass-loss events in LBV and SN progenitors 
are consistent. Note the 60 yr timescale is only an upper limit 
as we have not obtained a well-sampled temporal dataset of 
high-resolution spectra. Hence the variations are likely to be 
on shorter timescales still consistent with S Doradus variations. 
Even the change into a single absorption component (374 d) 
mimics the behaviour of LBVs such as AG Car which show sim- 
ilar variations on timescales of years (Fig. 3). 

For SN 2005gj the line profiles, wind velocities and wind pe- 
riodicity are remarkably reminiscent of LBVs and provide sig- 
nificant support for a direct link between supernovae and lumi- 
nous blue variables. This may have a huge impact on our under- 
standing of massive star evolution since current theory predicts 
that LBVs, which are in a short-term, intermediary phase of evo- 
lution, should not undergo core-collapse. 
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